Superconducting vortex avalanches, voltage bursts, and vortex plastic flow: efl'ect of 
the microscopic pinning landscape on the macroscopic properties 
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Using large scale simulations on parallel processors, we analyze in detail the dynamical behavior 
of superconducting vortices undergoing avalanches. In particular, we quantify the effect of the 
pinning landscape on the macroscopic properties of vortex avalanches and vortex plastic flow. These 
dynamical instabilities are triggered when the external magnetic field is increased slightly, and are 
thus driven by a flux gradient rather than by thermal effects. The flux proflles, composed of rigid 
flux lines that interact with 100 or more vortices, are maintained in the Bean critical state and 
do not decay away from it. By directly determining vortex positions during avalanches in the 
plastically moving lattice, we find that experimentally observable voltage bursts correspond to the 
pulsing movement of vortices along branched channels or winding chains in a manner reminiscent 
of lightning strikes. This kind of motion cannot be described by elastic theories. We relate the 
velocity field and cumulative patterns of vortex fiow channels with statistical quantities, such as 
distributions of avalanche sizes. Samples with a high density of strong pinning sites produce very 
broad avalanche distributions. Easy-flow vortex channels appear in samples with a low pinning 
density, and typical avalanche sizes emerge in an otherwise broad distribution of sizes. We observe 
a crossover from interstitial motion in narrow channels to pin-to-pin motion in broad channels as 
pin density is increased. 

PACS numbers: 64.60.Ht, 74.60.Ge 
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I. INTRODUCTION 

Dissipative extended systems that are slowly driven to 
marginally stable states produce complex and novel dy- 
namics. Systems characterized by avalanche dynamics, 
in which energy dissipation occurs in sudden bursts of 
collective activity, are of particular interest. Avalanche 
behavior has been studied in many systems, including 
granular assemblies magnetic domains |'^, charge 

density waves |^ , fluid flow down inclines [^|0 1 , and flux 
lines in type-II superconductors jl^-jl^]. 

Microscopic information can reveal why some regimes 
of material parameters produce broad distributions of 
avalanche sizes, while others do not. Such informa- 
tion is generally unavailable, however, and only limited 
macroscopic changes in system configurations can be ob- 
served. Further, experimentally tuning microscopic pa- 
rameters and recording detailed microscopic information 
about the dynamics is difficult. Numerical simulations 
allow exact control of microscopic parameters, such as 
pin strength and density, and provide both very precise 
microscopic information, such as individual vortex mo- 
tion, and macroscopic information, such as voltage sig- 
nals. Comparison of simulation and experiment has been 
hampered because the extreme numerical demands of ac- 
curate models has caused virtually all theoretical explo- 
rations of these systems (see, e.g., Ref. Q) to employ 
discrete dynamics governed by simple rules. With recent 
advances in parallel processing, however, it is now possi- 
ble to study much more realistic and complicated models 
using molecular dynamics (MD) simulations. In this pa- 



per, we present a continuous, MD simulation of supercon- 
ducting samples containing a critical state of very slowly 
driven vortex lines that undergo avalanching behavior. 
We find that the density of pinning sites plays an im- 
portant role in producing or suppressing broad distribu- 
tions of avalanches. When the pin density is low, favoring 
easy-flow vortex channels, characteristic avalanche sizes 
appear. For higher densities, no unique channels form, 
and distributions remain very broad. 

Flux penetrates a type-II superconductor in the form 
of discrete quantized vortices that repel each other and 
are attracted by defects in the superconducting material. 
A gradient in the vortex density develops and drives vor- 
tices into the material. A balancing pinning force holds 
the vortices in a metastable state, known as the criti- 
cal state or Bean state . As an external field is slowly 
increased, additional flux lines enter the sample and occa- 
sionally cause large disturbances. Although experiments 
can use local changes in flux to detect the motion of these 
vortices pl|jr^ , p^ -p^ , at present the vortex motion can- 
not be directly imaged over long enough time scales to 
permit a statistical characterization of the motion. Thus, 
computationally generated information on vortex move- 
ments is of great interest. 

To investigate dynamical instabilities producing cas- 
cades of flux lines in superconductors, we have performed 
extensive MD simulations on parallel multiprocessors us- 
ing a wide variety of relevant parameters that are diffi- 
cult to continuously tune experimentally, such as vortex 
density n.^, pinning density Up, and maximum pinning 
force (or strength) fp. We do not observe a parameter- 
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independent universal response to perturbations in our 
samples, but instead find a rich variety of behaviors in 
which all of these parameters play an important role. 
This ranges from the collective motion of vortex chains 
dominated by pin-to-pin transport to the appearance of 
very narrow interstitial channels, where vortices flow be- 
tween pinning sites, in agreement with recent Lorentz 
microscopy results [ p^ . Computer simulations 19-2l|] 
are a valuable tool for the analysis of the microscopic 
spatio-temporal dynamics of individual flux lines in su- 
perconductors and its relation with commonly measured 
macroscopic averages. With simulations, interactions be- 
tween a plastic vortex lattice and rigid pinning sites can 
be easily examined p^ ]. 

We use a T = MD algorithm to perform large-scale, 
detailed simulations of many superconducting samples. 
Each sample has one of three pinning strengths fp and 
one of three pinning densities rip. Five different combi- 
nations of fp and rtp, spanning a wide variety of possible 
pinning configurations, were considered in this work. We 
study dynamical^ rather than thermal, instabilities; thus, 
thermally-activated flux creep p^ ] does not occur in our 
system, and all avalanches are driven solely by the com- 
petition between the vortex gradient and pinning forces. 
Each avalanche is triggered by the addition of a single 
flux line to the system. After the avalanche ends and the 
system reaches mechanical equilibrium, another flux line 
is added to the system. That is, during each avalanche, 
no vortices are added to the sample. The work presented 
here is distinct from previous MD simulations in sev- 
eral ways. Our detailed study involves more than 10^ 
avalanches for each combination of fp and n^, recorded 
using an extremely large number of MD time steps (10* 
hours on an IBM SP parallel computer). This allows us 
to construct reliable, statistically significant distributions 
from the avalanches. Our more realistic two-dimensional 
model also employs a much longer vortex- vortex interac- 
tion range than previously used pO| . 

This paper is organized as follows. In section II, a de- 
tailed description of the numerical simulation is given. 
Section III contains our definition of a vortex avalanche. 
In section IV, we present a variety of images of the vor- 
tex avalanche events, highlighting common features such 
as vortex motion in winding chains and the pulse-like 
nature of the avalanches. In section V, statistical dis- 
tributions are analyzed. Features in these distributions 
are related quantitatively to the microscopic parameters 
and dynamics of the system. Sections VI and VII con- 
tain brief comparisons of our work to recent experiments. 
Finally, we summarize our results in Section VIII. 



II. SIMULATION 

We model a transverse two-dimensional slice (in the 
x-y plane) of an infinite zero-field-cooled superconduct- 
ing slab containing rigid vortices that are parallel to the 



sample edge (H = Hi). The sample is periodic in the 
y direction only, and there are no demagnetizing effects. 
An external field is modeled by the presence of flux lines 
in an unpinned region along one edge of the sample. This 
field is very slowly increased by adding a single vortex to 
the unpinned region each time the sample reaches a state 
of mechanical equilibrium. Vortices enter the supercon- 
ducting slab under the force of their mutual repulsion and 
pass through a pinned region 24A x 26A in size, where A is 
the penetration depth. They form a fiux gradient natu- 
rally due to their own interactions and give rise to the 
critical current dB/dx = 27r Jc/c. Experiments employ a 
flux profile which, on average, does not change with time 
inside the sample [|ll|. To model this, we remove vortices 
from the simulation when they exit the opposite end (i.e., 
the right edge in our figures) of the pinned region. 

The vortex- vortex repulsion is correctly represented by 
a modified Bessel function, Ki{r/X], cut off beyond r — 
6A where the force is quite small |Q. The vortices also 
interact with a large number of non-overlapping short- 
range attractive parabolic wells of radius = 0.15A. The 
maximum pinning force, /p, of each well in a given sample 
is uniformly distributed between f™^^/^ and f^^^. The 
pinning strength in our samples varies from very weak 
pins {fp < 0.3/o) to very strong pins {fp > 2.0/o), with 
maximum pinning values given by f™'^^ = 0.3/o, l.O/o, 
or 3.0/0. Here, 
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The samples studied here contain one of three differ- 
ent densities of pinning sites, rip, corresponding to three 



different numbers of pinning sites, Np-. 



5.93/A^ 



{Np = 3700), rip = 2.40/A2 {Np = 1500), orn^ = 0.96/ 
{Np = 600). Although the pinning potential is fixed in 
time, the energy landscape produced by the moving vor- 
tices evolves continuously. 

The overdamped equation of vortex motion is 
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where the total force fi on vortex i (due to other vortices 
f™, and pinning sites f"^) is given by 



fc=i 



A 



(2) 



Here, 9 is the Heaviside step function, is the location 
of the iih vortex, is the velocity of the ith vortex, r^?''' 
is the location of the fcth pinning site, is the radius of 
the pinning site, Np is the number of pinning sites, Ntj 



is the number of vortices, 



(r^ — r[f^)/|ri — r^*''' |, and we take 77 = 1. We measure all 
forces in units of /o = ^q/9i'k'^\^, magnetic fields in units 
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of $o/A^, and lengths in units of the penetration depth 
A. The number of vortices forming the Bean state varies 
from sample to sample, ranging from a minimum value of 
Ny w 240 for the sample with a low density of strong pins, 
Up — 0.96/A^ and f™'^^ = 3.0/o, to a maximum value of 
Ny « 1700 for the sample with a high density of strong 
pins, np = 5.93/A2 and f^^"" = 3.0/o. Throughout this 
paper we examine the regime 0.06 ^ riy/np ^ 0.82, where 
riy is the vortex density. 
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FIG. 1. Magnetic flux density profiles B{x) for each 
sample studied. The profiles were obtained using 
B{x) = {24\)-^ J^^^ dyB{x,y). The area < a: < 4A, 
< y < 24A, to the left of the dotted line, is the unpinned 
region through which vortices enter the sample. This mimics 
the external field; the dashed line indicates the field strength 
in this region. The presence of bulk pinning in the sample 
(located in the region 4A < x < 30A, < y < 24A) provides 
a barrier for flux entry and exit at the interface a: = 4A. The 
five proflles correspond to flve samples with three different 
densities of pinning sites, rip, and three different uniformly 
distributed pinning strengths, fp. Starting from the top 
(very strong pinning) to the bottom, we have filled triangles: 
/™ = 3.0/0, rip = 5.93/A^ Np = 3700, N^ fa 1700; open 
squares: f^'''' = 3.0/o, Up = 2.40/A^ Np = 1500, iV„ « 1000; 
open diamonds: /"^'^ = l.O/o, Up = 5.93/A^ A^p = 3700, 
iV„ « 700; asterisks: f^'''' = 3.0 fo, rip = 0.96/A^ A^p = 600, 
N^ ^ 500; and plus signs: j;^'"' = 0.3/o, rip = 5.93/A^ 
Np = 3700, iV„ ^ 240. Note that the slope of B{x), i.e., 
Jc{x), is somewhat larger towards the right edge of the sam- 
ple where the flux density is lower and the effective pinning 
is larger. The average slope is not altered by avalanches since 
the majority of the vortices in the sample do not move during 
an avalanche. 



III. VORTEX AVALANCHES 

The flux profile in the sample (see Fig. |l|) is maintained 
in a quasi-magneto static state, with only a single vortex 
added to the left side of the sample every time the sam- 
ple reaches mechanical equilibrium ||2^]. Although this 
resembles granular sandpile experiments that add a sin- 
gle grain at a time to the apex of the pile |jl|,||, there 
are very important differences between granular systems 
and our vortex system: the vortex interaction has a much 
longer range, incorporating the effects of up to 100 near- 
est neighbors; inertial effects are negligible in this highly 
overdamped system; and the disorder is quenched (i.e., 
pinning sites are fixed) and tunable (that is, the effective 
pinning strength can be varied with field). 

To gauge the amount of vortex motion occurring in a 
sample at each MD time step, we use the average vor- 
tex velocity Wav, given by the sum of the magnitude of 
the velocity of each vortex Vi divided by the number of 
vortices in the sample Ny, 
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1=1 



(3) 



We add a new vortex only when significant motion has 
stopped, as indicated when 



(4) 



where Wth is a low threshold value. The signal r;av is plot- 
ted in Fig. 1^ for a sample with a high density of strong 
pinning sites and a threshold level vth = 0.0006/o/ry. 
Each time a new vortex is added to the sample, a very 
small peak appears in Wav. Larger peaks are produced 
during the avalanches that are occasionally triggered by 
the addition of a vortex, with the largest peaks gen- 
erated by a combination of either a few vortices mov- 
ing rapidly or a large number of vortices moving slowly. 
Many avalanche events produce Wav signals that consist 
of several peaks clustered together, indicating that the 
avalanche disturbance does not propagate through the 
sample at constant v^-v, but moves in repeated pulsing 
waves 1 25 . This burst-like behavior is a result of the 
combined effects of the vortex density gradient and the 
two-dimensional nature of the vortex motion. 

Most avalanche disturbances start at the outer sample 
edge (left edge in the figures) where the vortex density 
is highest. When a vortex inside the sample is displaced 
from its pinning well due to a small increase in the ex- 
ternal field, pinned vortices ahead of it prevent it from 
moving down the gradient. This is in contrast to gran- 
ular sandpiles, where sand grains can move in the third 
dimension to avoid such blockages. A string of vortices is 
depinned in a domino- like effect, with each moving only 
far enough to depin the next vortex and then stopping. 
This pulsing motion continues until the forces on all vor- 
tices in the sample are once again below the threshold 
depinning force. 
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FIG. 2. Plot of the average vortex velocity t^av, shown in normal (a) and expanded (b) views. The box in (a) indicates the 
region expanded in (b). This quantity resembles voltage signals induced by moving vortices in actual experimental samples 
[11]. An additional vortex is added to the left edge of the sample only when Vav falls below the threshold (background) value 
of Vth = O.OOO6/0 /r). Every period of avalanche activity can be clearly resolved. Certain velocity bursts composed of several 
closely spaced peaks in Uav correspond to single avalanches in which several vortex chains move consecutively in a manner 
resembling lightning strikes. In these avalanches with multiple peaks in Uav, the disturbance does not propagate through the 
sample at constant Wav in one pulse, but moves in repeated pulsing waves. This signal is from a sample containing 1700 vortices 
interacting with a high density of strong pinning sites: Up = 5.93/A^, f^^^ = 3.0/o. 



By adding only one vortex at a time, we probe the 
system in the limit of a zero current of incoming vor- 
tices. Large values of the current do not produce well- 
defined individual avalanches, so it is important that a 
given avalanche die away conipletdy before; adding the 
next vortex. Adding many vortices simultaneously, or at 
a very fast rate, can significantly modify the state of the 
system. The information we learn about a rapidly driven 
system tells us not about the original state but about the 



(possibly drastically) modified state. Use of an effectively 
infinitesimal perturbation allows us to avoid altering the 
nature of the original system while still gaining informa- 
tion about it, something that is important in other sys- 
tems as wcill. For instance, in classical mechanics, a very 
small perturbation is applied to a particle to test the sta- 
bility of its orbit. In electromagnetic systems, a small test 
charge is added to systems in order to probe their elec- 
trical screening or dielectric properties. In many-body 
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physics, a very small electric or magnetic field is applied 
to systems to determine their diamagnetic response. 



A. Voltage bursts and intermittency 

The time-evolution of the velocity signal plotted in 
Fig. ^ resembles that observed for a variety of differ- 
ent physical systems in the so-called intermittent regime 
1^,^. In particular, the general features of the veloc- 
ity signal in Fig. |^ are similar to the fluid velocity at a 
point in the interior of the container in Rayleigh-Benard 
convection experiments ||2^ in the intermittent regime. 
In the vortex and fluid cases, the somewhat regular ve- 
locity signals are intermittently interrupted by velocity 
"bursts" of finite duration which occur at seemingly ran- 
dom times. This section explores the conceptual analo- 
gies and differences between the velocity bursts produced 
by vortices moving in superconductors and the velocity 
bursts observed in intermittency. A quantitative and de- 
tailed comparison is beyond the scope of this section, and 
will be presented elsewhere. 

Let us first summarize the intermittent transition 
to chaos jS. More information can be found in 
Refs. p6| , ^7[ . Intermittency refers to a signal that al- 
ternates in an apparently random manner between long 
regular phases, known as laminar phases or intermis- 
sions, and relatively short irregular bursts, called chaotic 
bursts. The frequency of these bursts increases with an 
external parameter that we call e here. Thus, intermit- 
tency provides a continuous route from regular, burst-free 
motion, for e = 0, to chaotic motion, for large enough e. 
For very small values of e, there are long stretches of 
time, called laminar phases, during which the dynamics 
is regular, remaining very close to the burst-free e = 
fixed point. 

Let us now consider a marginally stable Bean state. 
Take e to be the rate at which vortices are added to the 
system (i.e., the driving rate). If e = 0, the dynamical 
system of vortices is at a fixed point in both position and 
velocity spaces. In particular, its average velocity Uav is 
zero. In the language of nonlinear dynamics, the laminar 
phase now has an infinite duration and the motion is 
regular; that is, it is a fixed point of zero velocity. 

When vortices are added to the system at a very small 
rate e, such as one vortex added every 10^ MD steps, 
some vortices in the system rearrange their positions. 
In this case, the average vortex velocity remains very 
close to the v^v = laminar phase value, up to small 
oscillations produced by vortex rearrangements. In other 
words, if e > and very small, the velocity of the system 
remains quite near its Wav = fixed point (see, for in- 
stance. Fig. 50(b) of Ref. |^). This dynamics near the 
fixed point, or "inside the laminar region", does not con- 
tinue forever. Eventually, the system exhibits a burst; 
that is, Uav(0 suddenly increases to a larger value, pro- 
ducing a crest of apparently random shape, or a "chaotic 



burst". After a period of time that is typically short, 
this burst suddenly comes to an end, and the system is 
"reinjected into" the laminar region, exhibiting its usual 
dynamics near the fixed point. When vortices are added 
to the system at an increasingly large rate e, the average 
vortex velocity displays more frequent bursts away from 
the Wav = fixed-point laminar phase value. 

For chaotic systems in the intermittent regime, the 
length of time the system spends around its fixed point, 
and the length of time the system spends in the chaotic 
region, are both unpredictable. If the same computer run 
is repeated with initial conditions that vary slightly, the 
details of the dynamics will not be reproduced p6,p7[. 
These and many other results have been obtained for 
chaotic systems in the intermittent regime. Analogous 
calculations have not been made for vortices in the Bean 
state. It would be interesting to calculate, among other 
things, the average time the system spends around its 
fixed point (i.e., in the laminar region) as a function of 
the driving rate e. It is unclear at this point how these 
quantities compare for the vortex and standard nonlin- 
ear systems. To compute these quantities for vortices is 
far more complicated than for the usual systems. Fluid 
flow, modeled by the three Lorentz equations, and other 
intermittent systems (e.g., nonlinear RLC circuits) can 
be described effectively by few degrees of freedom. The 
problem of vortex avalanches requires solving very many 
degrees of freedom, which greatly complicates the calcu- 
lations. These will be considered elsewhere. 



IV. ANALYSIS OF THE SPATIOTEMPORAL 
EVOLUTION OF VORTEX AVALANCHES 

The spatiotemporal propagation of avalanches through 
a sample are of great interest. The analysis presented 
below provides insight into which microscopic properties 
of a sample are most important in determining the na- 
ture of the resulting avalanches. Few experimental stud- 
ies in any avalanche system are able to view individual 
avalanches. Experiments involving vortex avalanches so 
far have been able to resolve only a macroscopic result of 
the avalanche, such as flux leaving the sample, and have 
not been able to determine the path the vortices followed 
as they moved through the sample. Our simulations can 
easily image vortex avalanches, revealing how changing 
the microscopic pinning parameters affects the nature of 
the avalanches. 



A. Spatial configurations of vortex avalanches 

We first determine the path an avalanche follows 
through the sample by identifying the vortices involved 
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FIG. 3. Snapshots of vortices participating in an avalanche event. Hero, the initial positions of vortices that wore dcpinncd 
during the avalanche are represented by filled circles. Vortices that remained pinned are represented by crosses. Vortex motion 
is towards the right of the figure, and the entire 26A x 24A sample is shown. Panel (a) is from a sample with Up = 5.93/A^ 
and Z™"" = 3.0/0; panels (b), (c), and (d) arc from a sample with lower pinning density: rip = 2.40/A^ and f^^^ = 3.0/o. (a) 
shows a large event. The stationary vortices in and around the flow path arc more strongly pinned than the moving vortices 
and provide barriers to the flow, (b) is a typical chain-like event. Most events are this size or smaller, (c-d) show how a 
characteristic channel for motion can change slightly over time [e.g., from (c) to (d)]. Infrequent large events rearrange the 
vortices in (c) and alter the interstitial pinning caused by vortex-vortex interactions, resulting in the new channels shown in 
(d). In other words, the channels in (d) are produced by a vortex lattice rearrangement after the channels in (c) dominated the 
transport over several avalanches. Note the gradient in the vortex density, with a higher density on the left side of the figure. 
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in representative avalanche events. A vortex is consid- 
ered an avalanche participant if it is depinned, indicated 
when its displacement di during the time interval be- 
tween adding new vortices to the sample is greater than 
the pinning diameter 2^p, 

d^ > 2^p . (5) 

This distinction is necessary since most vortices in the 
sample are displaced very slightly during an avalanche, 
but only a small number depin and are transported down 
the gradient. Because our simulation operates in the 
regime n„ < Up, all vortices remain trapped at pinning 
sites when an avalanche is not occurring, except in the 
case of very low pin density when some vortices remain at 
interstitial sites created by the repulsion from surround- 
ing vortices. 

A single avalanche event is depicted in each panel of 
Fig. m The entire sample is shown in each case, with 
the initial positions of vortices that were depinned dur- 
ing the avalanche marked with filled circles. The vor- 
tices that remained pinned are indicated with crosses. 
Rare large events, such as that illustrated in Fig. ||(a) 
for a densely pinned sample, involve a broad region of 
the sample. These events occur only after a considerable 
amount of strain has accumulated in the vortex lattice. 
In the most frequently occurring events, winding chains 
of vortices move from pinning site to pinning site, with 
each vortex moving into the site vacated by a vortex to 
its right, as in Fig. ||(b). The chains extend down the 
flux gradient in the x-direction, and are not perfectly 
straight but wind in the ^-direction. The amount of 
winding increases as the pinning strength fp decreases. 
The chains do not appear at the same location during ev- 
ery avalanche, and chain size may vary: in some events a 
chain spans the sample, while in others the chain contains 
only three or four vortices. 

When the pinning density is high, moving vortex 
chains are equally likely to form anywhere in the sample. 
As the pinning density is lowered, however, the probabil- 
ity that the movement will occur in a certain well-defined 
channel at one location in the sample becomes very high. 
The position of this channel varies from sample to sam- 
ple, but all samples with low pinning density contain 
such a channel. An example of a "preferred channel" 
is illustrated in Fig. ||(c,d). Here the vortices follow the 
same winding path during several consecutive avalanches. 
These meandering vortex channels display branching be- 
havior or form small loops around stronger pinning sites. 
As a result of the low pinning density in the sample shown 
in Fig. ^(c,d), interstitial pinning is highly impoTtant. In- 
terstitial pinning occurs when a vortex is held in place 
only by the repulsion from surrounding core-pinned vor- 
tices. Since the strength of the interstitial pinning is 
significantly weaker than that of the surrounding core 
pinning sites in a sample, as in Fig. |^(c,d), the result- 
ing avalanches follow easy-flow paths composed almost 
entirely of interstitially pinned vortices. The occasional 



large events in such a sample rearrange the interstitial 
pinning landscape, leading to the formation of a new set 
of flow channels that will persist for a period of time until 
the next large event occurs. For example, the persistent 
flow channel of Fig. ^(c) was altered into the channel of 
Fig. 1(d) by a large event. 

A small packet of flux never actually moves from one 
end of the sample to the other in a single avalanche, re- 
gardless of whether vortices move in fixed channels or in 
constantly changing paths. Instead, movement is trans- 
mitted from vortex to vortex, with an individual vortex 
rarely moving more than one to two pinning sites away 
from its former location during an avalanche. Thus, the 
disturbance crosses the sample, but a vortex does not, 
and the time span of a typical avalanche, r, (measured 
below) is much shorter than the time required for a single 
vortex to traverse the sample. 

B. Snapshots of the velocity field 

The images analyzed above show that avalanche dis- 
turbances involve several possible types of moving bun- 
dles or chains of vortices, but give no information about 
the time scales for vortex motion. Thus, to examine the 
avalanche dynamics, we consider the velocities of indi- 
vidual vortices at several instants during an avalanche 
and present a series of snapshots of the vortex move- 
ments during a single avalanche. We find that avalanche 
disturbances propagate neither instantly nor at constant 
speed through the sample. 

Figure ^ shows an avalanche moving through a portion 
of a sample with a high density of strong pinning sites. 
An arrow of a length proportional to the instantaneous 
vortex velocity marks the initial position of each mov- 
ing vortex. The velocity of the vortices near the right 
edge of the sample is greater than that near the left edge 
due to the vortex density gradient. The uniform time 
interval separating each panel was chosen as th, the time 
required for a vortex to hop from one pinning site to 
the next, so that most vortices complete their motion 
within one panel of the flgure. The disturbance in this 
avalanche propagates through the sample, with most of 
the vortices involved stopping exactly one pinning site 
to the right of their initial positions. Thus, the major- 
ity of the motion occurs inside the sample and cannot 
be detected by probing only vortices exiting the sample. 
Those vortices that do exit have initial positions within 
one to two lattice units from the sample edge, and we 
do not observe vortices moving many pinning sites to the 
right during a single avalanche in order to leave the sam- 
ple. We therefore find that during the largest avalanches, 
when the greatest number of vortices exit the sample, all 
movement and flux exit occurs in a wide region of the 
sample. Fig. |^(a) shows an example of a wide avalanche. 

The event pictured in Fig. ^ has a medium- length life- 
time. Events with longer lifetimes often consist of more 



7 




FIG. 4. Consecutive snapshots of the vortex velocity field in a 25A x 5A region of a 26A x 24A sample with a high density of 
strong pinning sites, Up = 5.93/A^, /^^'^ = 3.0/o. Each moving vortex is indicated by an arrow whose length is scaled by the 
velocity of the vortex; the remaining vortices are indicated by small crosses. All vortices sit in pinning sites (not shown) when 
not in motion. A vortex is considered "moving" if it is dcpinncd. The remaining vortices arc not completely motionless, but 
shift very slightly inside the pinning sites. The disturbance propagates from the dense left edge of the sample to the relatively 
less dense right edge. The vortices in the rest of the sample (not shown) were not depinned. The time interval between 
snapshots is the typical time th required for a vortex to hop from one pinning site to another. The illustrated motion is typical 
of a medium-sized avalanche in this sample. 
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FIG. 5. Continuous lines indicate the paths vortices follow through a portion of the sample over an extended period of time 

covering many avalanches. Vortices enter through the left edge and move towards the right edge. Dots represent vortices, and 
opcu circles represent pinning sites. An 8A x 8A region of the 26A x 24A sample is shown, (a) rip = 0.96/A^, f^^^ = 3.0/o. (b) 
Up = 2.40/A^ /^^'^ = 3.0/0. (c) Up = 5.93/A^ /^^'^ = 3.0/o. (d) Up = 5.93/A^ /^^'^ = 0.3/o. Vortex motion in the sample 
with Up = 5.93/A^, Z™*" = l.O/o resembles (c) and (d). In (a) and (b), strong pinning occasionally causes segments of the vortex 
path to run towards the top and bottom of the figure, transverse to the flux gradient. The presence or absence of easy-flow 
channels is strongly dependent on pin density. The channels present in (a) and (b) lead to avalanches with characteristic sizes 
Na and lifetimes r. Samples with higher pinning density do not have isolated easy-flow channels and produce very broad 
distributions of avalanche sizes. Because of the higher pinning strength, avalanches in (c) have higher typical vortex velocities 
Vi and shorter typical event lifetimes r than those in (d), resulting in tighter, less wandering vortex paths. 



than one pulse of motion; that is, the average velocity avalanche. 
Va.v exhibits several spikes or oscillations during a single 
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C. Cumulative pattern of vortex flow channels 

The spatial configurations of the avalanches we observe 
are not strongly affected by the sample pinning parame- 
ters, with all samples producing vortex motion in wind- 
ing chains that propagate through the sample during one 
or more pulses. The pinning parameters are important, 
however, when we consider whether chains moving in con- 
secutive avalanches are concentrated in one area of the 
sample, or whether they are evenly distributed through- 
out the sample over time. By plotting the vortex trajec- 
tories with lines, and drawing these lines for an extended 
period of time covering many avalanches, we can identify 
the cumulative pattern of vortex flow channels for differ- 
ent pinning parameters. Heavily-traveled regions of the 
sample are easily distinguished by a concentration of tra- 
jectory lines. An examination of small regions of samples 
shows that the typical flow channel pattern varies with 
pinning density. 

As shown in Fig. |^(a-b), interstitial channels of easy 
flow for vortices develop in samples with a low den- 
sity of strong pinning sites. The flow of flux lines in 
Fig. ||(a) involves only the much more mobile interstitial 
vortices moving plastically around their strongly pinned 
neighbors, indicated by the fact that the vortex trajec- 
tories form narrow paths that never intersect the pin- 
ning sites. Similar behavior has been directly observed 
recently by Matsuda et al. where interstitial, chain- 
like avalanche flow was seen around a few strongly pinned 
vortices [Q. As the pin density is increased, pin-to-pin 
vortex motion becomes more important and the intersti- 
tial channels become less well defined [Fig ^(b)] until at 
the highest pin densities, Up = 5.93/A^, interstitial flow 
is no longer observed, and the vortices always move from 
pinning site to pinning site. This is illustrated in Fig. ^(c- 
d) for a high density of pins of two different strengths. 
Note the absence of distinct isolated flow channels. 



it is a distance dp/2 from each pinned vortex, and exerts 
a force on each equal to 

/mMpoint ^ -f^l ■ (6) 

For interstitial motion to occur, we require 

j-(p) _ fvv 
Jmin J midpoint " 

Thus, 

/mill ^ •^o (2a) 

For each of the pinning densities used in our simula- 
tion, we can determine the minimum pinning strength 
/mill that would permit interstitial motion. In samples 
with a high density of pinning sites, Up — 5.93/A^, as 
in Fig. |(c,d), we find f'^^l^ = 4.7/o. Since the actual 
pinning forces in these samples are fp ^ 3.0/o, we do 
not expect interstitial motion to occur, and the sim- 
ulations confirm that only pin-to-pin motion occurs in 
these samples. Samples with a lower pinning density of 
= 2.40/A2, as in Fig. |(b), have /^t = 2.8/o. In 
this case, with fp ^ 3.0/o, occasional interstitial motion 
may occur near the strongest pins, and a small amount of 
interstitial motion is observed in the corresponding sim- 
ulations. In samples with the lowest pinning density of 
Up = 0.96/A2, as in Fig. |(a), = I.6/0. This condi- 
tion is easily met by a large portion of the pins, which 
have fp ^ 3.0/0, a-nd exclusively interstitial vortex mo- 
tion is expected and observed in the simulations. 



D. Derivation of the criteria for the presence of 
interstitial vortex motion 

In the limit of strong pinning, we expect intersti- 
tial vortex motion to occur whenever a vortex moving 
between two pinned vortices exerts a maximum force 
/mMpoint sach pinned vortex that is less than the pin- 
ning force fp. This occurs at the midpoint between the 
two pinned vortices. If this condition is not met, one 
of the pinned vortices will be depinned, and the inter- 
stitial vortex might be trapped by the vacated pinning 
site. Given the pinning density Up of a sample, we can 

iv) 

calculate the minimum pinning strength f^(^ required to 
allow interstitial motion. Assume that two vortices are 
pinned in adjacent wells separated by the average dis- 

— 1/2 

tance dp — Up between pinning sites. If an interstitial 
vortex passes directly between these two pinned vortices. 



V. STATISTICAL DISTRIBUTIONS 
CHARACTERIZING VORTEX AVALANCHES 

The strikingly different behaviors of the vortex 
avalanches are quantified here by identifying the depen- 
dence of the distributions of several quantities on pinning 
parameters. These quantities include: t, the avalanche 
lifetime; Na, the number of vortices participating in the 
avalanche; Nf, the number of vortices exiting the sam- 
ple during an event; dtot, the total vortex displacement 
occurring in the avalanche; as well as di and Vi, the in- 
dividual vortex displacements and velocities during an 
avalanche. 



A. Avalanche lifetime 

A frequently used characterization of an avalanche is 
its total lifetime, r, or the interval of time during which 
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the avalanche occurs. In our simulation, r is equal to 
the interval between perturbations of the system by the 
addition of a new vortex. Recall that here, the Bean 
state is always driven in a quasimagnetostatic mode: a 
flux line is added, vortex positions in the sample shift 
and may produce an avalanche, and the next flux line 
is added only after the vortex lattice reaches mechanical 
equilibrium. We expect the avalanche lifetime to depend 
on both the pinning strength and the pinning density. 

The pinning strength fp influences avalanche lifetimes 
by determining the speed of vortex motion in the sample. 
We find a relationship between the pinning strength fp 
and the speed w,; of an individual vortex by considering 

the behavior of a vortex at position Fj immediately after 

(p) 

it is depinned from a well at position r j. • If we assume 
that the vortex barely has enough energy to escape the 
well, the total force on the vortex just before it depins is 
close to zero: 



r + fpr 







where 



/,^|fr(|r.,.H^p)| 



(8) 



(9) 



is the maximum pinning force at the edge of the parabolic 



well, Yik = Yik/\Yik\ = (r,; 



.(p)| 



and Tik is 



the distance between vortex i and the center of the fcth 
parabolic well. Thus, 



fp^ik • 



(10) 



The pinning force has an abrupt cutoff at the pinning 
radius ^p, so when the vortex moves off the pinning site, 
the vortex-pin force f"^ suddenly falls to zero while the 
long-range vortex-vortex force t™ changes only negligi- 
bly. The resulting force on the vortex is 



Thus, the vortex velocity is 

|f^l ^ fp 



(11) 



(12) 



where Vc is a characteristic velocity associated with the 
pinning strength. Vortices in samples with stronger pin- 
ning move faster when depinned than vortices in samples 
with weaker pinning. 

In order to directly relate vortex velocities to avalanche 
lifetimes, we consider the distance an individual vortex 
moves during an avalanche event. A vortex normally 
hops from one pinning site to an adjacent site during the 
event. The distance it travels is simply dp, the average 
distance between pinning sites, 



1 



(13) 



If we define a "string" avalanche as an event during which 
each vortex in a chain extending the length of the sample 



hops from one site to the next, we can estimate a "string" 
lifetime. First, we designate the time that a vortex with 
just enough energy to depin spends moving between pins 
as the hopping time, th- 



th 



J] 



fp 



(14) 



G^-s p52^, so using 



For example, in PbIn, rj ^ 3.3 x 10 
fp — 3.0/0, n-p = 5.93/A^, and A ~ 65 nm gives th ^ 15 
fs. Next, we estimate Nh, the number of vortices that 
must hop in order for the avalanche to span the sample: 



(15) 



where is the sample length. If we assume that the 
vortices hop one at a time in sequence, we obtain an 
estimate of the "string" avalanche lifetime. 



thNh 



fpy/n^ 



(16) 



For instance, using the value of th for PbIn given above 
and taking the vortex density to be Uy = 2.5/A^ gives a 
lifetime of r^st ~ 0.6 ps. The relationship between Tost 
and fp and Up is not simple because the vortex density 
rijj = ny(H, fp,np). For a given field strength, however, 
it is clear that avalanche lifetimes decrease as pinning 
strength or density is increased. This is confirmed in the 
plot of the distribution P(t) in Fig. ^(a-b), where we see 
that samples with strong dense pinning have significantly 
shorter avalanche lifetimes than samples with weaker pin- 
ning [Fig. ^(a)] or samples with lower pinning density in 
which interstitial pinning is important [Fig. |^(b)]. All 
histograms have been smoothed as described in Rcf . [Q . 

In order to determine what fraction of the avalanches 
have lifetimes on the order of the estimated "string" life- 
time Test, we scale the avalanche lifetimes by Tost, and 
plot P^r/Tost) in Fig. ^(c-d). For the dense strong pin- 
ning case of Fig. ^c), most avalanche lifetimes are shorter 
than the estimated lifetime Tost {t /Test < 1) as a result 
of two factors. First, several vortices hop simultaneously, 
as in Fig. ^, rather than hopping one after the other, as 
assumed in our estimate. Second, although Fig. || illus- 
trated several system-spanning avalanches, many events 
do not involve enough vortices to span the entire sam- 
ple but contain only a short chain of vortices moving a 
short distance. The majority of avalanches in samples 
with strong pinning, J™^'^ = 3.0/o, contain less than ten 
vortices traveling in short chains that originate in regions 
of high vortex density but do not extend into areas with 
lower vortex density. 
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FIG. 6. Distributions of avalanche lifetimes, r. The left- 
most panels (a,c,e) correspond to our highest density of pins, 
Up = S.QS/A'^, and differing pinning strengths: filled trian- 
gles, Z™'"' = 3.0/0; open diamonds, Z™"" = l.O/o; plus signs, 
jmax _ Q 3j|^ The right panels (b,d,f) correspond to sam- 
ples with our strongest pins, J™'"' = 3.Q/o, and differing pin- 
ning densities: filled triangles, rip = 5.93/A^; open squares, 
np = 2.40/A'^; asterisks, Up — Q.96/A^. (c-d) Distributions 
scaled by the estimated lifetime. Test. Note the appearance of 
characteristic avalanche times indicated by arrows in (d) at 
t/tcsi ~ 0.8 and t/tcsi ~ 1.3. (e-f) Distributions scaled by 
the hopping time, th- In (e) the hopping time between pin- 
ning sites is th =6, 18, and 57 MD steps, respectively. Also 
in (e), notice the power law behavior P{T/th) ~ 
for small r/th over two decades for the sample with a high 
density of strong pins (filled triangles). In (f) the hopping 
time between pinning sites is th =6, 17, and 53 MD steps, 
respectively. The arrow in (f) at r/th ~ 30 indicates a local 
maxima in P{T/th) for the two samples with lower density of 
pins. Panel (g) schematically indicates the notation used in 
the figures of this work. 
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FIG. 7. Distributions of the number of vortices, Na, mov- 
ing during an avalanche event. The left panels (a,c,e) cor- 
respond to samples with a high pin density, Up = 5.93/A^, 
and difi'ering pinning strengths: filled triangles, Z™^" = 3.0/o; 
open diamonds, f^"-"" = l.O/o; plus signs, f^"-"" = 0.3/o. The 
right panels (b,d,f) correspond to samples with strong pins, 
jmax _ 3 ojg^ and differing pinning densities: filled trian- 
gles. Up = 5.93/ open squares, Up = 2.40/A^; asterisks. 
Up = 0.96/A^. (c-d) Distributions scaled by the average num- 
ber of vortices, Nh, that must hop for the avalanche to cross 
the sample. In (d), note the appearance of a characteristic 
value in the latter two samples at Na/Nh ~ 1. (e-f) Dis- 
tributions scaled by the average number, A'^„, of vortices in 
the sample. For samples with a large density of strong pins, 
the leftmost curve in (e) indicates that most avalanches in- 
volve only a tiny fraction of the vortices, resulting in plastic 
transport that occurs in brief, choppy bursts. 

As the pinning force is lowered, the chains of mov- 
ing vortices become longer on average, leading to an in- 
crease in the average value of r/Tost- Only samples with 
a reduced pinning density, shown in Fig. ^(d), produce 
frequent "string" avalanches. In these samples, nearly 
all avalanches occur in a single easy-flow channel that 
spans the sample, as in Fig. @(a), or in two main chan- 
nels as in Fig. ||(b). As a result, we observe an increased 
likelihood for the "string" avalanche in the form of an 
increase in the distribution function near r/rost ~ 1- 
This increase is most pronounced for the lowest pin den- 
sity, Up — 0.96/A^ [Fig. 11(d)], where the channel is most 
strongly established. Since the channels in such sam- 
ples wind significantly in crossing the sample, the num- 
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ber of vortices moving during the avalanche is greater 
than the estimated value Nh, and the observed increase 
in P(T/Test) for Up = 0.96/A^ falls at r/rcst ~ 1-3, rather 
than at r/rcst ~ 1- 

Figures ^(e-f) show the scaling of the avalanche life- 
times by the hopping time t^. In this case, the distribu- 
tions Pjr/th) for high pinning density {up = 5.93/A^), 
in Fig. |6|(e), collapse and can be approximated for small 
T by the form 




might be expected to have small values of Na- It is, how- 
ever, the width of the avalanche disturbance rather than 
the magnitude of the flux gradient that is important in 
determining Na- In samples with strong pinning, ener- 
getic avalanches occur, but avalanche width is suppressed 
since strongly pinned vortices on either side of the mov- 
ing chain are not depinned during the short interval when 
vortices are moving very rapidly down the steep vortex 
density gradient. As the pinning weakens, avalanches 
become wider when weakly pinned vortices adjacent to a 
slowly moving chain depin and join the motion down the 
shallow gradient, resulting in more events with relatively 
large values of Na- 



7 - 1.4 . 

Distributions generated by samples with lower pinning 
densities, in Fig. |6|(f), do not scale and cannot be repre- 
sented by the same form. This indicates the importance 
of the nature of avalanche propagation. For samples with 
a high density of pins, pin-to-pin vortex motion occurs 
throughout the sample, with all regions of the sample 
participating in an avalanche at some time. Once the 
pinning density decreases, however, interstitially pinned 
vortices appear and dominate avalanche transport. The 
vortices always follow the same paths, visiting as many 
interstitial sites as possible when moving through the 
sample, and introducing a characteristic avalanche life- 
time. 



B. Number of vortices in each avalanche 

Avalanches can be characterized according to the num- 
ber of vortices displaced during the event. This quantity 
can be directly observed through computer simulations, 
but at present must be inferred from experiments [ pT| . 
We define the number of vortices, Na, that were actively 
moving participants in each avalanche as the number of 
vortices that were depinned. Figure 0(a-b) presents the 
corresponding distribution P(Na). For high pinning den- 
sity, as in Fig. 0(a), we can approximate the distribution 
for small Na as 

PiNa) ^ N-P . 

We find that (3 decreases as the pinning strength de- 
creases: 



/3 



f^-^ = 3.0/0 
/p-- - l.O/o 
= 0.3/o 



Smaller values of (3 indicate a relative increase in the 
frequency of large avalanches compared to small ones. 
The trend in (3 might appear counterintuitive since the 
smaller fiux gradient present in systems with lower pin- 
ning strength fp produces less energetic avalanches that 
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FIG. 8. (a-b) Distributions of the total distance, dtot, 
moved by all vortices in the sample scaled by the width of 
the sample, L^- (c-d) Distributions of the average distance, 
dav, traveled by a vortex during an avalanche event scaled 
by the pin diameter, 2^p. The left panels (a,c) refer to sam- 
ples with high pin density, — 5.93/A^, and differing pin- 
ning strengths: filled triangles, f^^^ = 3.0/o; open diamonds, 
f^'''' = l.O/o; plus signs, /^'^'^ = Q.3/o. The right panels (b,d) 
refer to samples with strong pinning, /^^'^ = 3.0/o , and dif- 
fering pinning densities: filled triangles, Up = 5.93/A^; open 
squares, Up — 2.40/A^; asterisks, Up — 0.96/A^. 

Vortices fiow only in certain channels when Up is low, 
leading to the appearance of a characteristic value in 
P{Na) [arrow in Fig. 0(b)] and making it impossible to 
describe P{Na) by a power law form over the entire range 
of observed Na values. The formation of avalanches with 
large values of Na is impeded in such cases because all 
available interstitially pinned vortices are already par- 
ticipating in the channel, and the neighboring strongly 
pinned vortices cannot be depinned. 

The arrow in Fig. ^(b) indicates a value of Na cor- 
responding to the number of vortices in the channels 
-^channel- Sincc the channels span the sample, we ex- 
pect -/Vchannoi to bc closc to the uumbcr of vortices Nh 
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that must hop to cross the sample in a straight chain, 
where — Lx\/n^. To verify this, we plot the scaled 
quantity P{Na/Nh) in Fig. @(c-d). We find in Fig. 0(d) 
that samples with low pin density rip, in which inter- 
stitial pinning is important, have A'channoi ~ Nh, as is 
expected if all moving vortices form a single chain across 
the sample. Specifically, increases in P{Na/Nh) fall at 
Na/Nh w 1 for Up = 2.40/A2 and Na/Nh ~ 1.1 for 
Up — 0.96/A^. The winding of the channels in samples 
with Up = 0.96A^ causes the value of Na corresponding 
to channel motion to be slightly larger than N^. For 
samples with high pinning density. Fig. ^c) indicates 
that the most linear part of P{Na/Nh) is produced for 
Na ^ Nh, or for events that did not cross the entire sam- 
ple. Broad avalanches are more than one vortex chain 
wide and have Na/Nh > 3. These avalanches are more 
common in samples with smaller /p, as indicated by the 
rightmost curve in Fig. ^c), where samples with weak 
pinning provide the largest avalanches for a fixed proba- 
bility density P{Na/Nh). 

To determine how effectively avalanches transport vor- 
tices across the sample, we scale Na by the total number 
of vortices Ny present in the sample, and plot the re- 
sult in Fig. ^e-f). Only a tiny fraction of the vortices 
participate in a typical event when pinning is strongand 
dense, as indicated by the leftmost curve in Fig. 0(e). 
The avalanches in this sample are also very short lived, 
as already shown in Fig. ||(a), so this sample is best char- 
acterized by plastic transport that occurs in brief, choppy 
bursts. On the other hand, for a high density of weak pin- 
ning sites, we find that it is possible for a significant por- 
tion of the vortices in the sample to collectively move in 
an avalanche, as in the rightmost curve of Fig. |(e). Vor- 
tex motion in these long-lived avalanches is less plastic, 
with small adjacent portions of the vortex lattice gradu- 
ally sliding forward at different times. The samples with 
low pin density, shown in Fig. 0(f), have their transport 
dominated by the single easy-flow channel that develops. 
We find that the fraction of vortices contained in this 
channel increases as pinning density decreases since the 
interstitial channel winds to a greater degree in the sam- 
ple with fewer pins and involves a correspondingly larger 
number of vortices. 



C. Total vortex displacement 

Another measure of avalanche size, uniquely available 
through simulation, is the total displacement dtot of all 
vortices in the sample during the event. 



dtot — di ) 



(17) 



length with dtot > L^;, we plot P{dtot/Lx) in Fig. |(a- 
b) . Occasionally we observe avalanches with the surpris- 
ingly large dtat/Lx ^10. These large values reflect the 
cumulative effect of vortices that experience very small 
displacements inside the parabolic pinning sites rather 
than moving to a new pinning site. For instance, if 800 
vortices each move a distance di ~ ■Cp/lO = 0.015A, the 
total displacement recorded would be 12A. Such very 
small vortex displacements play a very important role in 
transmitting stress throughout the vortex lattice. Fig- 
ure ^(a) indicates that P{dtot/ L^) is very robust against 
variations in pinning strength. Excluding large values of 
dtot/Lx, the distributions can be approximated by the 
form 



P 



dtot 



We find that 5 decreases as the pinning strength de- 
creases: 

1.7, /--'^ = 3.0/0 
1.2, /-- = 1.0/o 
1.1, /-- = 0.3/0 

Figure ||(b) indicates that P{dtot/ L^) is weakly depen- 
dent on rip. The effect of the channels appears as an en- 
hancement of the distribution near dtot/Lx = 1, marked 
by an arrow in Fig. H(b). 

We can compare the relative amount of vortex trans- 
port in different samples by considering the average dis- 
tance dav moved by a vortex during an avalanche. 



dgrv — 



N, 



dtot 



(18) 



Since most vortices do not become depinned, we plot 
Pida.v/'2£,p), where 2^p is the diameter of the pinning trap. 
The plot of P{ds,v/2^p) in Fig. ||(c) shows that as the pin- 
ning strength fp is reduced, dav/2^p increases, indicating 
that individual vortices are likely to travel further during 
an avalanche. This occurs both because a large fraction 
of vortices in samples with weak pinning participate in 
the avalanche, as demonstrated in Fig. ^(e), and because 
the average spacing between vortices increases since vor- 
tex density decreases with decreasing pinning strength. 
In Fig. |8|(d) , we find that reducing the pinning density rip 
can either increase or decrease d^^. As Up decreases, the 
strength of the interstitial pinning sites decreases, allow- 
ing individual vortices to travel further. Not as many vor- 
tices are moving, however, since the motion is restricted 
to well defined channels. These two competing effects 
enhance P(dav) at small values of dav, but reduce it at 
large values. 



where di is the displacement of each vortex. In order to 
find the probability of an avalanche that spans the sample 
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FIG. 9. (a-b) Distributions of the distance, di, traveled 
by individual vortices during all the avalanche events. Char- 
acteristic displacements are marked with arrows, (c-d) 
Distributions scaled by the inter-pin distance, dp. Only vor- 
tices that moved more than l.Odp were directly participating 
in the avalanche. The rest were shifting very slightly in pin- 
ning wells and transmitting stress through the lattice, (e-f) 
Distributions P{di/dp) for only those vortices that traveled a 
distance di < dp, i.e., vortices that were not depinned. (g-h) 
Distributions P{di/dp) for only those vortices that traveled a 
distance di > dp, i.e., vortices that were depinned. The left 
panels (a,c,e,g) correspond to samples with a high density 
of pins. Up = 5.93/ \^ , and differing pinning strengths: solid 
Une, Z™^" = 3.0/0; dashed line, = l.O/o; heavy solid 

line, f^'^^ — 0.3/0. The right panels (b,d,f,h) correspond 
to samples with strong pinning, /"^^ — 3.0/o, and differing 
pinning densities: solid line, rip = 5.93/A^; dot-dashed line. 
Up = 2.40/ dotted line. Hp = 0.96/ 



D. Individual vortex displacements 

The microscopic information available in our simula- 
tion permits us to calculate individual vortex displace- 
ments. In Fig. ^ we plot P{di), where di is the distance 



a single vortex is displaced in an avalanche. In every dis- 
tribution, we find that vortex motion can be described in 
terms of a characteristic size dt, indicated by arrows in 
the figure. This size roughly corresponds to the average 
distance between pinning sites dp, 



dh 



1 



(19) 



and appears because vortices participating in an 
avalanche typically hop from one pinning site to an ad- 
jacent pinning site. Thus, in Fig. |^(a), d/j falls at the 
same value for three samples with equal pinning den- 
sity Tip, while in Fig. ^(b), dh increases as Up decreases. 
The vortices in the sample can thus be divided into 
three categories: those that did not leave a pinning site, 
di < 2^p < dh, those that hopped from one pinning site 
to another, di ~ dh, and those that hopped more than 
one pinning site, di > dh- 

We first consider those vortices that hopped one pin- 
ning site. If we scale the distributions by the inter-pin 
distance dp, as in Fig. ^(c-d), we find that the character- 
istic distances in P{di/dp) fall at di/dp ~ 1 (marked with 
an arrow), with a second smaller increase at di/dp ~ 2 
occurring in samples where vortices occasionally hopped 
two pinning sites. 

Next we examine the motion of vortices shifting very 
slightly inside pinning wells. In Fig. |(e-f), we focus on 
vortices that did not leave a pinning site, di < dp, and 
present P{di/dp). We find that, for small values of di/dp, 
these distributions can be described by the form 



where: 



P 



1.4 for 



di 



5.93/A2 



P ■ 



1.4 for rip = 5.93/ A^ and f^"" 
-- 2.40/A2 and fp^^ = 3.0/o; p 



and /■» 



and /--'^ = 3.0/0; 
= 1.0/0; P ~ 1-4 for 
1.2 for Up = 0.96/\^ 
3.0/0; and p w 0.9 for Up = 5.93/X^ and 
jmax _ 0.3/0. The two samples with smaller values of p 
also have the lowest critical currents Jc, as determined 
from the slope of the flux profiles in Fig. l]. 

A simple two-dimensional argument shows that p is 
expected to be similar (and close to one) for all sam- 
ples, since only vortices that experience extremely small 
displacements in pinning sites affect the value of p. If 
we add an additional vortex to an already existing ar- 
rangement of vortices, the repulsive force, /extra, from 
this vortex will cause each nearby vortex to move a very 
small distance Ui in its pinning well in order to reach 
a new equilibrium position. Neglecting higher order ef- 
fects, we find that the very small displacement Ui of the 
ith vortex is determined by the distance from the added 
vortex and by the form of the force /extra- We take the 
distance from the perturbing vortex to be r, and assume 
that the additional force is applied during a time St. The 
vortex-vortex interaction force gives 
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/e: 



(20) 



If we consider the limit of small r, r ^ A, we can use the 
asymptotic form of the Bessel function to write 



/e: 



H A 



87r2A3 r 
We have (for < < A) 

/extra 



5u{r) 



V 



-St 



6t A 
T] Stt^A^ r 



(21) 



(22) 



To obtain a distribution, wc find the number of vortices 
N{r) located a radius r from the perturbing vortex. Due 
to the cylindrical symmetry, this can be written simply 
as 



SN{r) = 2'Krny Sr , 



(23) 



where Sr is an infinitesimal displacement. We want an 
expression for the slope of our distributions, 



d 



■InSN . 



d{\n{Su)) 
Thus, we write: 

\nSN{r) = ln{2nnvSr) + Inr 

d\nSN{r) _ 1 
dr r 



\n.Su{r) « In 



dlnSu 



dr 



_3^\ 

1 

r 



Inr 



Therefore, we find 
dln^iV 



din Su 



= -1 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



For small displacements {Su < <C A), this predicts 
p ~ 1, in general agreement with our observed values 
(p« 0.9- 1.4). 

The argument presented above is independent of the 
shape of the pinning potential, and assumes that the vor- 
tex density n„ is constant throughout the sample. The 
presence of the Bean critical state; makes this assumption 
inaccurate. It can, however, be shown that including the 
critical state does not significantly affect the expected 
value for the slope. Wc describe the field in the sample 
according to the Bean model, 



B = He- 



27rX 



(30) 



where Hcxt is the externally applied field. The field can 
also be written as 



Thus, 



n„(r) 



B = n„$o • 



gext - 2-KJcr/c 



(31) 



(32) 



The number of vortices located at a distance r becomes 

SN{r) = 2,, (-^ext-27r Jer/c) ^^ ^33^ 
$0 

We then have 



\nSN{r) = In 



2'iTSr 



In r + ln(ifext - 2-KJcr/c) (34) 



d\TiSN{r) _ 1 -2-kJc/c 
dr r He^t -2Tr Jcr/c 



(35) 



The external field i?ext can be written in terms of the 
critical current Jc and the distance the field has pene- 
trated into the sample, r*: 



2TrJrr* 



We obtain 
dlnSN 



din Su 



1 — {r*/r) 



(36) 



(37) 



We have r* » r, since the field has penetrated the entire 
sample length ^ A. We therefore write 



dln^iV 
din Su 



-1 



1 



1 - (r/r*) 



-1 



r 



(38) 



Thus, including the Bean gradient introduces only a neg- 
ligible correcting term r/r* to the expected slope, since 

r/r* < 1. 

It is interesting to consider what slope value would be 
obtained if the vortex- vortex interaction were that found 

in a thin film, rather than the form for a thick slab of 
material used above. In this case, the additional force 
from a vortex is 



/e: 



1 



giving 



Su(r) — ^St 



In Su(r) ~ In ^ 2 In r 

V 



(39) 



(40) 



(41) 
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dlnSu 2 
dr r 



din SN 
din 5u 



(42) 



(43) 



Thus, the distribution of individual displacements in a 
thin film should differ from that produced by a slab of 
material {p ^ 0.5 versus p ~ 1). 

In Fig. ^(g-h), we focus on the displacements of vor- 
tices that were actively participating in avalanches, and 
plot P{di/dp) for vortices with di > dp. We find that de- 
creasing the pinning strength leads to a very slight overall 
increase in the distance traveled by a moving vortex, as 
indicated by the slightly higher likelihoods for moving 
larger distances in Fig. ^(g). This is a result of the fact 
that, in general, avalanche disturbances are larger and 
longer lived in these samples, allowing individual vor- 
tices the opportunity to travel two or more pinning sites 
in an avalanche. Decreasing the pinning density, as in 
Fig. |9|(h), also affects the likelihood that a vortex will 
move a certain distance. 
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FIG. 10. (a-b) Distributions of the number Nf of vortices 
falling off the edge of the sample, (c-d) Distributions scaled 
by the average total number Af„ of vortices in the sample. The 
left panels (a,c) correspond to samples with a high density of 
pins, rip = 5.93/ and differing pinning strengths: filled 
triangles, Z™"" = 3.0/o; open diamonds, J™'"' = l.O/o; plus 
signs, f^^^ = 0.3/0. When Nf is not very large, the slope of 
the curve in each case is roughly —2.4. The right panels (b,d) 
correspond to samples with strong pinning, J™^" — 3.0/o, and 
diff'ering pinning densities: triangles, Up = 5.93/A^; squares. 
Up = 2.40/A^; asterisks, Up — 0.96/A^. For small values of 
Nf, the slope of the curves in (b) is roughly —2.4, —3.4, and 
—4.4, respectively. 



E. Vortices leaving the sample 

The number of vortices Nf that exit the sample during 
an event can be directly compared with values obtained 
in experiments, both in vortex systems ||llHl3[ and in 
sandpiles ||l|-||. In Fig. 0(a-b) we plot PjjVjf. Events 
with Nf = are not shown. The relatively small size of 
our sample resulted in a smaller data set for this quantity 
than for all other quantities considered in this paper. If 
we approximate the distributions for small Nf by the 
form 

PiNf) ^ NJ" , 
we find that all samples with high pinning density have 
a - 2.4 . 

As the pinning density Up is decreased for constant /p, a 
increases, with a ~ 3.4 for Up — 2.40/A^ and a ^ 4.4 for 
Up — 0.96/A^. A larger a indicates a relatively smaller 
likelihood for the occurrence of events with large Nf. 

When we examine the region from which vortices exit 
the sample, we find that flux can exit from any location 
in samples with high pin density, as seen from the rela- 
tively uniform coverage of trails along the right (outer) 
edge of Fig. ||(c) and Fig. ||(d). Even after one location 
has been depleted by a large avalanche, other areas along 
the sample edge still contain enough vortices to remain 
active in large and small events while the depleted regions 
refill. As the pinning density is decreased and interstitial 
channels develop, avalanche paths become highly con- 
strained and flux exits in only a few locations that have 
difficulty building up enough stress to permit events with 
large Nf to occur. For example, with the reduced num- 
ber of flux paths in Fig. ||(b), events with large Nf are 
rare, and thus a is large: a ^ 3.4. Fig. ||(a) illustrates 
a small region of the extreme case of a sample with only 
one exit channel, for which Up = 0.96/A^ and a has the 
larger value a ~ 4.4. Large values of Nf occur only when 
this single channel is filled by vortices under high stress. 
Since all smaller events happen through this same chan- 
nel, however, stress in the channel is relieved by small 
events before it can build to a high level. Thus, the 
likelihood that the channel will produce a very large Nf 
is extremely small. As expected, the cut-off of P{Nf) 
for this sample falls at a much lower Nf than that of the 
other samples, as seen in the leftmost curve in Fig. p^(b). 

Strong pinning effectively confines the avalanche dis- 
turbance to a narrow channel and prevents movement 
from spreading throughout the sample. Weak pinning 
permits collective motion in larger regions of the sam- 
ple. These effects are more pronounced for samples with 
a lower density of pinning sites. This is highlighted by 
a comparison of how effectively individual avalanches re- 
move vortices from different samples. In Fig. p^(c-d) , 
we scale Nf by the number of vortices in the sample N^. 
Most vortices in a sample with strong pinning cannot 
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exit; those that do exit move in isolated chains. As the 
pinning strength is weakened, however, a greater percent- 
age of the vortices are able to exit in events that consist 
of the collective motion of several adjacent chains. This 
is illustrated by the rightmost curves in Fig. [ic|(c). 



F. Individual vortex velocities 

Due to the presence of a non-uniform magnetic flux 
gradient that is not of the ideal Bean form, the velocity 
of an average vortex is a function of position (see Fig. ||). 
Therefore, to construct a distribution of individual vor- 
tex velocities Vi that is not affected by the gradient, we 
select a narrow region of the sample where the gradient 
is essentially constant and find P{vi) for vortices in this 
region. The individual velocities Vi are then scaled by 
the characteristic velocity Vc, given by 



th 



(44) 



where dp is the average distance between pinning sites, 
and tfi is the average hopping time, th = V^-p/ fp- Fig" 
ure |l^ presents the resulting distributions P{vi/vc)- For 
a high density of pinning [Fig. |ll|(a)], we find that pin- 
ning strength changes Vc but does not affect the form of 
P{vi/vc)- Note that vortices moving at speeds less than 
Vc do not have sufficient energy to escape from a pinning 
site and are merely experiencing very small displacements 
inside their respective pinning wells. The similar form of 
P{vi/vc) is then explained by the fact that all of the wells 
are parabolic. For the lowest velocities, we can approxi- 
mate P{vi/vc) with 



P 



where 



Different densities of pinning wells alter the form of 
P{vi/vc), as can be seen from Fig. |l](b). As the pin- 
ning density decreases, some vortices experience small 
displacements in interstitial wells while others remain in 
parabolic traps. The combined effect changes the overall 
distribution. 

If, instead of considering the distribution of individual 
velocities, we construct the power spectrum of the ve- 
locity signal Wav, we find that all of the samples with a 
high pinning density produce broad power spectra S{f) 
of the form S{f) ~ f'", where 1.5 < < 2.0. Samples 
with lower pinning density do not produce spectra of this 
form. The noise power spectra will be discussed further 
elsewhere lls^l. 
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FIG. 11. Distributions of the individual velocities, Vi, of 
vortices in a thin strip of the sample, scaled by the charac- 
teristic velocity Vc- The left panel (a) corresponds to samples 
with a high pinning density, Up = 5.93/ , and differing pin- 
ning strengths: filled triangles, f^^^ = 3.0/o; open diamonds, 
= l.O/o; plus signs, f^'"''' = 0.3/o. The right panel (b) 



corresponds to samples with strong pinning, /™ 
and difi'ering pinning densities: triangles, Up 
squares, Up — 2.40/A^; asterisks, Up — 0.96/A^. 



= 3.0/0, 
5.93/A^ 



VI. RECENT EXPERIMENTS ON VORTEX 
PLASTIC FLOW AND AVALANCHES 

In several experiments, the nature of the vortex move- 
ment has been inferred indirectly by using, for instance, 
voltage signals, resulting in several plausible scenarios of 
vortex motion including cylindrical bundles |17|| , clump- 
like bundles [ ^2|j3^ , or elongated bundles [t36[ . Due to 
the presence of a field gradient in the sample, we find 
that all collectively moving groups of vortices are ar- 
ranged in one or more long, narrow moving chains rather 
than clumped or cylindrical bundles. The existence of 
quasi-onc-dimensional vortex channels is suggested in 
Rcfs. 32,3^], where vortex motion is assumed to occur 
in straight paths that cross the sample. The chain-like 
channels that we observe resemble these suggestions, ex- 
cept that in our simulation the vortex paths often wind 
through the sample rather than following the shortest 
path down the gradient. In addition, our simulation in- 
dicates that the velocity of the moving vortices during a 
single avalanche is not constant but is pulse-likc, that is, 
with oscillations in the average velocity during the 
vortex avalanche. 

In the experiment by Field et al. pT[ |, flux lines ex- 
iting a superconducting tube were detected with a coil. 
Our average velocity signals, plotted in Fig. ||, strongly 
resemble the resulting voltage signals in Fig. 1 of Ref. 
|]ll| . Both simulation and experiment find distributions 
of avalanche sizes Nj that can be approximated by power 
laws with a range of exponents. The experiment in Rcf. 
||ll| finds a 1.4 to 2.2, while we find a ~ 2.4 to 4.4 using 
a variety of different samples. The behavior of the exper- 
imental sample therefore resembles our simulations with 
a very high density of pinning sites, rip = 5.93/A^, and no 
characteristic channels. This is reasonable since the den- 
sity of pinning in the experimental sample is very high. 
The pinning centers are grain boundaries with spacings 
randomly distributed between about 30 and 50 nm, while 
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Quantity 


rmax 
Jp 


= 0.3/0 


rmax 


= 1-0./0 


rmax 
Jp 


= 3.0/0 


Ref. [11] 


Ref. [13] 




7 ' 


■V. 1 /I 


7 ' 


•V. 1 /I 


7 ' 








Na 


/3 


- 0.9 


P 


- 1.0 


P 


- 1.4 






dtot / Lx 


S r 


- 1.1 




- 1.2 




- 1.7 






di/dp 


P ' 


- 0.9 


P ' 


- 1.2 


P ' 


- 1.4 






Vi/Vc 


a ' 


- 2.0 


a ' 


- 2.0 


a ' 


- 2.0 






Nf 


a 


- 2.4 


a 


- 2.4 


a 


- 2.4 


a ~ 1.4 to 2.2 


a ~ 1.7 to 2.2 



TABLE I. Best-fit slopes of the most linear region of distributions of several quantities. The form of the distributions in each 
case were: P(r/th) ~ (T/th)-\ P(iV,) ~ iV"", P(dtot/L.) ~ (dtot/i.)"', P(d,/dp) ~ (d./dp)-^ Piv^jv,) - {v^jv^y , 
P(^Nf) ~ ■ Notice that a multiplicative factor in the independent variable does not alter the slope. Thus, r and r/th 
produce distributions with the same slope. The samples listed here had a high density Up — 5.93/A'^ of pinning sites. Samples 
with /™'"' — 3.0/0 and lower densities rip of pins did not produce power law distributions. 



interstitial motion; from Eq. |7|; this occurs when 

1 



\2\{T)^. 



Writing t = T/Tc, 

A(T) 



A(0) \/r^' 



(45) 



(46) 



and the strength of the pinning force fp required for in- 
terstitial motion to occur increases with temperature: 



A is on the order of 400 nm, leading to a pinning density of expected to happen whenever the pinning force exceeds 
Up ~ 100/ A2. Thus, this sample has a higher pin density the minimum pinning strength f^}^ that would permit 
than our most densely pinned case, and it is therefore 
reasonable to expect the experimental exponents to be 
less than 2.4 based on the general trend of our results. 

Zieve et al. |l^] observed avalanches at higher field 
strengths, where the pinning forces are expected to be re- 
duced, and interpret their results to mean that the high- 
field regime produces large-scale, plastic rearrangements 
of the vortices, whereas the low-field regime generates 
elastic rearrangements. It is possible that there were un- 
detectable steps in the hysteresis loops they obtained in 
the low- field regime since the smallest flux jumps resolved 
by their apparatus were on the order of 0.1 kOe. If their 
sample contained either characteristic channels for flux 
flow or high strength pinning, only a small fraction of 
the flux inside the sample would be transported through 
narrow exit regions at these low fields. As the field den- 
sity increased, their sample would cross into the regime 
of weak dense pinning. In this shown by our sim- 

ulations, the sample can effectively transport much more 
flux, producing larger flux steps that can be resolved by 
the experiment. 

Finally, we consider very recent work by Nowak et al. 
[ p^ , who obtain different types of avalanche distributions 
as a function of temperature. At higher temperatures, 
they find a broad distribution for the amount of flux ex- 
iting the sample in each event. 



(47) 



P{Nf) ^ NJ" , 

with a ranging from 1.7 to 2.2, in agreement with both 
Ref. [0 and our results. At the lowest temperatures, 
characteristic avalanche sizes appear in the form of events 
involving a large number of vortices that may be system- 
spanning. These changes have been discussed in terms 
of thermal instabilities in the material [T^ , but they can 
also be considered from the standpoint of the presence 
or absence of channels for vortex motion. As we have 
seen in our simulations, channels are most likely to form 
when the vortices are able to move interstitially. This is 



At the lowest temperatures, /roq is small, and interstitial 
motion is possible in the sample. Channels of flux flow 
form, and a characteristic avalanche size appears in the 
distribution of avalanche sizes. As the temperature in- 
creases, /loq increases until some of the pinning sites in 
the sample are no longer strong enough to permit inter- 
stitial flow. In this case, pin-to-pin vortex motion will 
occur evenly throughout the sample, and the distribu- 
tion of avalanche sizes will broaden. We therefore expect 
a transition from pin-to-pin motion at higher temper- 
atures to interstitial channel flow at low temperatures, 
with a corresponding transition from broad distributions 
of Nf (for high T) to a characteristic value of N f (for low 
T). The transition with temperature in the nature of the 
distributions is experimentally observed in Ref. [[T^ . 

Table I presents a summary of our results for the ex- 
ponents of the power laws found in our avalanche distri- 
butions. Only the distribution of flux lines falling off the 
edge of the sample was measured experimentally [|ll|,0 . 
The other quantities are more difficult to measure ex- 
perimentally. However, novel flux imaging techniques 
(such as magneto-optical imaging, arrays of Hall probes 
|0,B8|, scanning Hall probes, and especially, Lorentz 
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microscopy) could make it possible to obtain some of 
these distributions experimentally. Indeed, flux-gradient 
driven vortex rivers similar to the ones described in this 
work have already been imaged using Lorentz microscopy 




VII. COMPARISON WITH OTHER AVALANCHE 
STUDIES 

Superconductors represent only one of the many sys- 
tems exhibiting avalanche behavior that have recently 
been studied. In this section, we briefly compare our 
work with a small sample of studies on avalanches in 
other systems, focusing on dynamical instabilities in dis- 
sipative extended systems which are very slowly driven 
towards (and not away from) marginally stable states. 
The literature on this subject is vast, and it is not the 
goal of this section to review it. 

The work presented in this paper differs from previous 
studies in several ways. First, most theoretical studies on 
avalanches (e.g., Ref. 10, |39j, and references therein) 
employ simple dynamical rules. Our system evolves ac- 
cording to realistic equations of motion and uses a re- 
alistic range of physical parameters. Second, instead of 
using ad hoc dynamical rules acting on discrete space 
and time variables, our simulations employ a dynamics 
which is continuous in both space and time. Third, most 
discrete dynamical rules proposed so far involve inter- 
actions among nearest-neighbor cells, or at most next- 
nearest neighbor cells. Every one of our vortices can in- 
teract with over 100 neighboring vortices, making possi- 
ble truly cooperative "cascades" in the marginally stable 
Bean critical state. This type of highly correlated motion 
is difficult to realistically model with discrete maps. 

Although most theoretical studies of avalanches have 
revealed no internal structure in the dynamics of the 
avalanches, Ref. describes a new series of ex- 

tremal (i.e., uniformly driven) models that produce 
avalanches with an interesting hierarchical structure of 
subavalanches within avalanches, not found in the earlier 
simpler models. These subavalanches appear to be simi- 
lar to the ones observed in our vortex avalanches in the 
form of velocity oscillations, as shown in Fig. ^. 

Avalanches exhibiting similar broad distributions of 
sizes and lifetimes have been observed experimentally 
in a variety of systems that are otherwise very differ- 
ent. In these systems, the movable objects can be vor- 
tices, grains, electrons, or water droplets. These interact 
through different types of forces, from very short range 
(hard-core) interactions for granular assemblies to longer 
range forces for electrons and vortices. The movable 
objects are driven by a variety of different forces (e.g., 
flux-gradient, Lorentz force, electrical current, gravity) 
and dissipate energy while they are driven (e.g., due 
to particle-particle collisions). These systems can ex- 
hibit both static and dynamic friction. For instance, the 



pinned state of vortices is the analog of static friction for 
grains, while the dissipative flux-flow regime is the ana- 
log of the dynamic friction seen in vortices. The inertial 
effects are very important in some cases, such as granular 
motion, and negligible in other cases, such as vortex mo- 
tion. The disorder in the sample can be frozen in, as in 
the quenched disorder produced by pinning sites in a su- 
perconductor, or dynamically evolving, as in avalanching 
granular assemblies. In spite of all these differences, each 
of these systems exhibits avalanches with broad distri- 
butions of sizes and lifetimes. Some of the distributions 
follow a power law over a limited range of values, re- 
sulting in similar exponents (as seen in Table II). Most 
non-superconducting systems produce exponents in the 
range from 2 to 2.5. Superconducting vortex avalanches 
can produce exponents that are significantly below 2. In- 
deed, the large variety of pinning landscapes present in 
superconducting samples allows the possibility of observ- 
ing a range of exponents from 1.5 to 2.2. 

Several words of caution are needed when comparing 
the exponents presented in Table II. First, the ranges over 
which power laws have been observed in these systems 
are small, typically covering between one and three or- 
ders of magnitude in the independent variable {x axis). It 
would be ideal to probe the response of each system over 
many more decades, but this is difficult to achieve exper- 
imentally. Even with realistic numerical simulations it is 
difficult to study a very large number of particles over 
a long period of time in order to probe large avalanche 
sizes and lifetimes. Second, the quantities displayed in 
Table II are not identical. In some cases, the measured 
"avalanche size" is the number of particles (e.g., grains, 
vortices, droplets) falling off the edge of the sample, while 
in other cases it refers to the actual number of particles 
that moved in the sample without leaving it. Third, the 
manner and rate at which particles are added to the 



Ref. 


Method 


Power law exponent 


[9] 


Water droplets 


1.9 


[7] 


Magnetic garnet films 


2.5 


[1] 


Granular pile 


2.5 


[42] 


Granular pile 


2.5 


[2] 


Granular pile 


2.2 


[6] 


Rice grains 


2.0 


[5] 


Ramped granular slides 


2.1 


[11] 


Superconducting vortices 


1.4 - 2.2 


[13] 


Superconducting vortices 


1.7 - 2.2 



TABLE II. Exponents of power laws observed in vari- 
ous experiments. [9] considered avalanches in a continuous 
medium. [7] observed magnetic domains. [1] worked with 
AI2O3 particles and beach sand. [42] used 3 mm iron or glass 
beads as well as 1 - 2 mm plastic beads. [2] placed 0.4 - 
0.8 mm Si02 sand in piles of varying sizes. [6] monitored a 
quasi-one-dimensional pile of rice. [5] imaged the surface be- 
havior of sand avalanches. [11] and [13] studied avalanches in 
"hard" (i.e., with srong pinning sites) superconductors. 
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system can affect the values of the power law exponents 
or even alter the functional form of the avalanche distri- 
butions. In some cases, incoming particles are added in 
the bulk, while in others, particles are added only at one 
edge. For instance, experiments in which incoming par- 
ticles are randomly sprinkled on the sample are driven 
in the bulk, while in experiments such as Ref. and 
PI , the particles are added on one edge only, resulting 
in different exponents. Simulations also produce differ- 
ent exponents for bulk driving and for driving on 
one edge 1^,1^. Fourth, comparisons among systems 
displaying avalanche dynamics are difficult to make be- 
cause of the significant differences among these complex 
systems, as noted above. 

Many questions still remain open, including why so 
many different systems exhibit broad distributions of 
avalanche sizes and lifetimes, what the origin of the power 
law behavior is, and why the values of the power law ex- 
ponents range between 1.5 and 2.5. In spite of consider- 
able efforts, a complete and convincing answer to these 
and other related questions is still lacking. 

VIII. SUMMARY 

We have quantitatively shown how the micro- 
scopic pinning parameters determine the nature of the 
avalanches that occur in superconducting samples driven 
very slowly by an increasing external magnetic field. By 
using large-scale MD simulations to monitor the vortices 
participating in avalanches, we observe motion along 
winding paths through the sample, and find that each 
vortex moves only one to two pinning sites during an 
avalanche rather than crossing the entire sample. Most 
avalanches are small and are contained inside the sam- 
ple. Thus, they cannot be detected with experiments 
that probe only vortices exiting the sample. 

Avalanches produce two kinds of vortex motion: pin to 
pin, and extremely small displacements inside the small 
parabolic wells. Only vortices that move from pin to 
pin directly participate in the avalanches, suddenly re- 
leasing accumulated stress in the vortex lattice through 
a succession of choppy bursts. The vortices that do not 
directly participate in a given avalanche experience ex- 
tremely small displacements (typically di <^ dp). These 
small shifts help to slowly build up and transmit stress 
throughout the vortex lattice. 

Pinning strength affects the average lifetime of a given 
avalanche by determining under how much stress the lat- 
tice is held. A high density of strong pinning holds the 
vortex lattice under a great deal of accumulated stress, lo- 
calizes avalanche disturbances, and is effective at retain- 
ing vortices inside the sample. The resulting avalanches 
are suppressed in width, involve rapidly moving vortices, 
and have very short lifetimes. Thus, avalanches in sam- 
ples with strong dense pinning are best characterized by 
plastic transport that occurs in brief, choppy bursts along 



narrow vortex paths Weak or sparse pinning pro- 

duces little build up of stress in the vortex lattice. The 
avalanche events in such samples involve vortices moving 
slowly in long-lasting events along much broader vortex 
channels. Higher vortex mobility in samples with weak 
pinning leads to larger total vortex displacement dav [p5| . 

The presence or absence of distinct channels for flow 
leads to a crossover from broad distributions to charac- 
teristic avalanche sizes. At low pin densities, avalanches 
pulse through the sample in narrow heavily trafficked 
winding channels composed of interstitially pinned vor- 
tices. As pin density increases, pin-to-pin vortex mo- 
tion dominates, the isolated channels disappear, and 
avalanches display a broad distribution over more than 
a decade. An important, and non-intuitive, result of our 
simulations is that the size of the critical current is not 
a good indicator of broad distributions. Indeed, samples 
with very different critical currents Jc (e.g., a sample with 
/ma'^ = 3.0/0, Tip = 0.96/A2, and - 0.03$oc/27r and 
a second sample with f^^^ — 0.3/o, Up — 5.93/A^, and 
Jc ~ 0.01<I>oc/27r) may have very different breadths to 
their distributions; the presence of unique channels, not 
the critical current, is the important factor. 

We note that avalanches in the Bean state cannot be 
characterized by universal avalanche distributions valid 
for all values of the pinning parameters. Although the 
Bean state is always critical, it does not always display 
avalanches with a lack of characteristic scale, and in those 
cases where avalanche distributions take the form of a 
power law, the exponent varies with pinning pUJT^ . Our 
results are consistent with experiments and explain the 
sample- and regime-dependence of recent experiments 
011 • 
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